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DETERMINING BOTH SOUND SPEED AND 
INTERNAL SOURCE IN THERMO- AND 
PHOTO-ACOUSTIC TOMOGRAPHY 


HONGYU LIU AND GUNTHER UHLMANN 

Abstract. This paper concerns thermoacoustic tomography and 
photoacoustic tomography, two couple-physics imaging modalities 
that attempt to combine the high resolution of ultrasound and 
the high contrast capabilities of electromagnetic waves. We give 
sufficient conditions to recover both the sound speed of the medium 
being probed and the source. 


1. Introduction 

In this paper we consider thermoacoustic tomography (TAT) and 
photoacoustic tomography (PAT), two coupled-physics modalities that 
attempt to combine the high resolution of ultrasound and the high 
contrast capabilities of electromagnetic waves. In the latter a laser 
pulse is used to probe a medium and in the former a lower frequency 
wave is used. As a result the medium is heated and produces an elastic 
expansion that in turns produces a sound wave that is measured outside 
the medium. This is used to hud optical properties of the medium in 
the case of PAT or electromagnetic properties in the case of TAT. These 
two imaging modalities and other coupled-physics imaging techniques 
have the potential of signihcant clinical and biological applications [3l 

Eiiniiig. 

We next present the mathematical formulation of the PAT and TAT 
problem described above. Let C be a compact set such that 
is connected. Let c{x) G L°°(M^) be such that supp{l — c) C fl, 
and c(x) > Cq for x E fl, where Cq is a positive constant. Consider 
the following Cauchy problem of the wave equation for u{x, t), (x, t) G 
X M+, 

f d‘fu(x,t) — Au(x,t) = 0 in X M+, 

< c\x) (1.1) 

\ u{x,0) = f{x), dtu{x,0) = 0 for x G 

where /(x) G L°°(M^) with supp{f) C fl. c(x) represents the wave 
speed in the body and /(x) represents an internal source produced 
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by the tissue volume expansion. The wave measurements are given by 

Af^cix,t) = u{x,t), (x, t) G X M_|_. (1.2) 

The inverse problem is to recover (/, c) by knowledge of Af^c{x,t). In 
what follows, for simplicity, (/, c) shall be referred to as a TAT/PAT 
conhguration. 

The uniqueness issue for the TAT/PAT problem has gained extensive 
attention in the literature. However, most of the available results in¬ 
tend to recover either / or c, by assuming that the other one is known. 
For the recovery of / by assuming that c is known, we refer to [1] and 
the references therein for the constant sound speed case; and [9] for the 
variable sound speed case. In [7j, the recovery of a smooth sound speed 
c by assuming that the source / is known was considered. Uniqueness 
was established in [7] with one measurement under the geometric as¬ 
sumption that the domain is foliated by strictly convex hypersurfaces 
with respect to the Riemannian metric g = where dx^ denotes 

the Euclidean metric. In [10], the instability of the linearized problem 
in recovering both the internal source and the sound speed was estab¬ 
lished. In [B] , stability estimate was derived in recovering / if there is a 
modeling error of the variable sound speed. In [2] , an interesting con¬ 
nection was observed between the TAT/PAT problem and the interior 
transmission eigenvalue problem. 

In this work, we show that in certain practically interesting scenarios, 
one can recover both the sound speed and the internal source by a single 
measurement. Our argument relies on the temporal Fourier transform, 
converting the TAT/PAT problem into the frequency domain. Then by 
using the low frequency asymptotics, we derive two internal identities 
involving the unknowns. The main uniqueness that we could derive 
from the two integral identities is that if the sound speed is constant and 
the internal source function is independent of one variable, then both of 
them can be recovered. We also point out some other uniqueness results 
in our study. Moreover, our argument requires the least regularity 
assumptions on the parameters as well as the domain. The rest of the 
paper is organized as follows. In Section 2, we present the main results 
of the current study. Sections 3 is devoted to the proofs. 

2. Statement of the main results 

We hrst introduce the admissible class of TAT/PAT conhgurations 
(/, c) for our study. Let (/, c), H and u{x,t) be given in fll.l|l . Our 
argument relies essentially on the temporal Fourier transform of the 
function u{x,t) defined by 

1 

u{x,k) = — / M(x,t)e*^* dt, {x,k) G x M+. (2.1) 

27 ! Jo 
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A TAT/PAT configuration (/, c) is said to be admissible if fl2.ip defines 
an distribution for k G (0, eo) with a certain eo G IP+. It is 

remarked that a sufficient condition to guarantee the admissibility of a 
TAT/PAT configuration (/, c) is that c is a non-trapping sound speed; 
see and the references therein for more relevant discussion. 

In what follows, we let 

g*fc|x| 

—TT for 7^ 0. 

A'k\x\ 

<h is the fundamental solution to —A — For the subsequent use, we 
also set 

for kl 7^ 0. (2.2) 

47r|x| 

Let Bp> := i?i?(0) denote a central ball of radius R G M+ containing hi. 
Moreover, for a function -0 G supported in Bji, we denote 

:= [ ^o{x-y)-i/jiy) dy, a; G (2.3) 

JR3 

We shall prove that 

Theorem 2 . 1 . Let {fj,Cj), j = 1,2, be two admissibleTAT/PAT con¬ 
figurations such that 

= Af^^c 2 {x,t), {x,t) G (912 X M+. (2.4) 

Then one has 

[ {cf‘^fi-cf‘^f 2 ){x)-v{x) dx = 0, (2.5) 

JR3 

where v(x) is an arbitrary harmonic function, and 


[ (1 - Ci^)(2/) • ^ \A^fi)iy) ■ ^o{x -y) dy 

9k3 

+ ^ / {cf^fi){y) -Ix-yl dy 
= / (1 - c^^)(2/) ■ ^"^(C 2 V 2 )(l/) ■ $o(a; - y) dy 

9r3 

+ ^ / {cf^f 2 ){y) -Ix-yl dy, 

OTT 

which holds for any x G dBji. 


( 2 . 6 ) 


It is noted that in fl2.5p and fl2.6p the integrations can be taken over 
12 or Br instead of since both (/i, ci) and (/ 2 , C 2 ) are supported in 
12. In the sequel, we let x = {xj)^^^ G Using Theorem 12.11 one can 
show that 
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Theorem 2.2. Let {fj,Cj), j = 1,2, be two admissibleTAT/PAT con¬ 
figurations. Then one has that 

cf‘^{x)fi{x) = cf‘^{x)f 2 {x) for a.e. X e Vt, (2.7) 

if either of the following two conditions is satisfied 

i). (p{x) := cf‘^{x)fi{x) — cf‘^{x)f 2 {x) is a harmonic function in 

Q; 

ii). ip{x) from item i) is independent of one variable, say Xj, 1 < 

J<3. 

Remark 2.1. By Theorem 12.21 one readily sees that if an admissible 
TAT/PAT conhgnration (/, c) is snch that either c~^f is a harmonic 
fnnction in hi, or c~^f is a fnnction independent of one variable Xj, 
1 < j < 3, then c“^/ is nniqnely determined by for {x,t) G 

(912 xR+. Using this result, one can further infer that if the sound speed 
c is known, then the internal source / is uniquely determined; whereas 
if / is known in advance, then the sound speed c\supp{f) is uniquely 
determined. 

Next, we present a practical scenario in uniquely determining both 
the sound speed and the internal source. 

Theorem 2.3. Let (/, c) be an admissible TAT/PAT configuration sup¬ 
ported in 12 such that c is constant, and f is independent of one variable 
Xj, 1 < j <?)■ Furthermore, we assume that 

f{x) > 0 for a.e. x G 12 and / f{x) dx > 0. (2.8) 

Jn 

Then both f and c are uniguely determined by Af^c{x,t) for {x,t) G 
(912 X M+. 

Remark 2.2. Assumption 12 .8 1 on the internal source / means that there 
exists an open portion of 12 on which / is strictly positive. Since / is 
generated by the absorbed energy, it is a practically reasonable condi¬ 
tion. 

By generalizing Theorem 12.31 we can further prove a bit more general 
unique determination result as follows. 

Theorem 2.4. Let (/, c) be an admissible TAT/PAT configuration 
such that 

c~^ix) = cf^{x) + 7“^Xs, a; G 12, (2.9) 

where Cb is a positive background sound speed which is supposed to be 
known in advance, and 'y is a positive constant denoting an anomalous 
inclusion supported in S C 12. Furthermore, we assume that both Cb 
and f are independent of the variable Xj, 1 < j < 3 and 

J A“^(c“^/)(x) • v{x) dx d ^ 


(2.10) 
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for a harmonic function v. Then both f and 7 are uniquely determined 
by Af^c{x,t) for {x,t) G dil x M+. 

From our subsequent argument, one can see that if Cb and S in 
Theorem 12.41 are taken to be 0 and and moreover v = 1, then 
Theorem 12.31 is a particular case of Theorem 12.41 


3. Proofs 


It is straightforward to show that u{x, k) satishes 
Au{x,k) + ^^u{x,k) = - — ^f{x), (x,/c) e X (0,eo). (3.1) 

C^[X) ZTT C^[X) 

Moreover, in order to ensure the well-posedness of the transformed 
equation fl3.ip . one imposes the classical Sommerfeld radiation condi¬ 
tion 

lim \x\(^^^^^^-iku{x,k)] =0. (3.2) 

|x|^+cx) \ d\x\ J 

Under the temporal Fourier transform, the TAT/PAT measurement 
data become 

A(x, k) = u{x, k)\onx{o,eo)- (3-3) 

The following lemma shall be needed. 


Lemma 3.1. Let u{x,k) G be the solution to fl3.1^ - 03.2p . 

Then u{x, k) is uniquely given by the following integral equation 


u{x, k) 


-k'^ f [1 - c ‘^{y)] 
./R3 


■u{y,k) ■$(x-r/) dy 


ik 


C ‘‘(y) ■ f(y) ■ Hx - v) dy, X 6 R’. 


(3.4) 


Moreover, as k ^ -|-0, we have 

Tik f 

u(x,k) = -— c~^{y)-f{y)-^o{x-y)dy + 0{k‘^), x e Br. (3.5) 


Proof. By the elliptic regularity (see, e.g, 0 ), we known that u G 
The integral equation fl3.4p is of Lippmann-Schwinger type, 
and it is uniquely solvable for u G C{Br) (cf. [1]). By using the fact 
that 

<h(x) = <Fo(a^) + C>(/c) as k ^ -fO, (3.6) 

uniformly for x G Br, one can verify fl3.5p by straightforward asymp¬ 
totic estimates. 

The proof is complete. □ 


We are in a position to present the proof of Theorem 12.11 
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Proof of Theorem \2.1{ Let Ui{x,k) and U 2 {x,k) denote the acoustic 
wave fields, respectively, corresponding to (/i, Ci) and (/ 2 , C 2 ). By fl2.4p 
and 03.31) . one clearly has that Ui{x,k) = U 2 {x,k) for {x,k) G dfl x 
(0, eo). Since both iti and U 2 satisfy the same equation (A + k‘^)u = 0 
in one readily has that 

Ui{x, k) = U 2 {x, k) for (x, fc) G x (0, eo), (3.7) 

and therefore 


Ul{x, k)\dBjix{0,eo) — U2{X, k)\dBiiX{0,eo) ■ 

Hence, by 03.41) and 03.Sp . we see that for x G OBb, 


(3,8) 


-k 

i 

= -k I 

i 

~ 2ti 


'Br 


[1 - Cl ^( 2 /)] • ui{y, k) ■ $(a; - y) dy 


Cl (y) ■ fi{y) -Hx-y) dy 


’Br 


'Br 


[1 - C 2 ^( 2 /)] • * 2 ( 2 /, k) ■ $(a: - y) dy 


(3.9) 


C 2 ( 2 /) ■ f 2 {y) ■ <h(a: - y) dy. 


'Br 


Next, by Lemma [3.11 one has 
ik 


Uj{x,k) = I Cj^{y)-fj{y)-%{x-y) dy+0{P), j = 1,2, (3.10) 
^^Jbr 

for k G M+ sufficiently small. Plugging O3.10p into 03.9p and using 
the series expansion for of <h(a: — 2 /), one can further show by 

straightforward calculations (though a bit tedious) that as fc —)■ +0 

k 


7 ^/ {y) ■ fi{y) ■ - y) dy 


27r 


'Rii 

i 

^ ' 

i 

IOtt^ 


Stt^ 


Cl {y)-fi{y)dy 


'Br, 


fcM (1 - cr^( 2 /)) ■ A-^(crVi)(2/) ■ $o(x - y) dy 

JBr 

/ ci'^iy)-fi{y)-\x-y\dy + 0{k^) 


'Br 


= A- / C 2 ( 2 /) • f 2 {y) • $o(a; - y) dy 


271 


'Br 


k 

Stt^ 


C 2 ( 2 /) ■ f 2 {y) dy 


'Br 


- — ■k^ 
271 

i 


'Br 


(I-C2 ( 2 /))-A (C 2 f 2 )iy) ■ ^o{x - y) dy 


IOtt^ 


k / C 2 (y)-f 2 {y)-\x-y\dy + 0{k'^), 


'Br 


which holds for any x G OBb- 


(3.11) 
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From fl3.1ip . one immediately has fl2.6p and for any x G OBr, 


'B. 


Cl {y)-fi{y)-%{x-y) dy = / C 2 {y)-f 2 iy)-^o{x-y)dy. (3.12) 


'Bf, 


We proceed to the proof of the orthogonality relation fl2.5p . and hrst 
note the following addition formula for <l>o(x —?/) with |a:| > \y\ (cf. pQ), 


47 r|x — y\ 


00 a 

a=0 l3=—a. 


2 a + 1 




(3.13) 


where x = x/\x\ and y = y/\y\, and Y^{x) denotes the spherical har¬ 
monics of order a G M U {0}, for {3 = —a,... ,a. By inserting fl3.13p 
into fl3.12p . together with straightforward calculations, one then has 


00 a 

YY 

a=0 ^=—OL 

00 a. 

xz 

0=0 I3=—ol 


1 


2 a -|- 1 
1 1 


• / q ‘^{y) ■ fi{y) ■ li/l" • Yi {y) dy 

JBr 


2 a -|- 1 i?“+^ 


■y!( 


X) 


'B, 


C2^(y) ■ f2{y) ■ iy, 

(3.14) 


which holds for any x G dBn. Since {F)f (•)}o=o,i,...,/ 3 =-a,...,o is a com¬ 
plete orthonormal basis of where denotes the unit sphere, 

one immediately has from 03.141) that 


'Bf, 


'Bf 


Cl (1/) ■/i(i/) • |i/|“ • 


C 2 (y)- f 2 {y)-\y\‘^ -Y^iy) dy, 


(3.15) 


for a = 0,1, 2,..., and {3 = —a, ..., a. We note that |r/|“y)f (^) for 
a = 0, 1 , 2 , ..., and (3 = —a, ..., a yield all the homogeneous harmonic 
polynomials. Hence, 03.15P readily implies 02.5p . 

The proof is complete. □ 


Proof of Theorem, \2.31 For case i), if <p(x) is a harmonic function in fl, 
one can take v{x) to be a harmonic function such that v{x) = <p(x) for 
x G n. Then by 02.5p in Theorem 12.11 one immediately has that (p = 0. 

For case ii), without loss of generality, we assume that y){x) is inde¬ 
pendent of X 3 for X = {xi,X 2 ,x^) G Let 

v{x) = X G M^, (3.16) 

where 

C = e + e = (6,6, 0) G = (0, 0, ^3^) G (3.17) 

satisfying ^1 + ^2 = is straightforwardly verihed that v{x) in 

03.16p - 03.17p dehnes a harmonic function. By taking v{x) in 02.5p to 
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be the one defined in fl3.16p . we have 


0 = / Lp{xi, X 2 )e^‘^'^ dx = / ip{xi,X 2 )e^‘’'^ dx 


>Br 



i^{x') dx' ■ 


/ 

-'{x3;(x',X3)G-Bh} 


dxs, 


(3.18) 


where x' = {xi,X 2 ) G and = {^ 1 ,^ 2 ) G In fl3.18l) . we have 
made use of the fact that ip is supported in Br. Clearly, from fl3.18p . 
one readily has 


dx' = 0, 


(3.19) 


which holds for any G M^. The LHS of 03.1911 is the Fourier transform 
of the function 9?, and hence one easily infers from 03.1911 that 9? = 0. 
The proof is complete. □ 


Proof of Theorem 1^.31 Without loss of generality, we assume that / is 
independent of the variable x^ and write it as f{x') for x' = (xi,X2) G 
By contradiction, we assume that there exists another TAT/PAT 
configuration (/(x'),c) supported in with c being a constant such 
that 

Af^c{x,t) = Aj-~{x,t) for {x,t) G dil x M+. (3.20) 

Since c“^/ — c“^/ is independent of the variable X3, we have by 
Theorem 12.21 that _ 

c-V = c-V- (3.21) 

Therefore, 



x-y\ dy 


Jbr 


x-y\dy 


(3.22) 


and ^ 

A-\c-^f) = A-\p-^f) (3.23) 

Next, by 02.6p in Theorem 12.11 together with the use of 03.2211 and 
03.23p . there clearly holds 


(c-2-h-^)/ A-\c-^f)iy)-%{x-y)dy = 0, (3.24) 

Jbr 

where x G BBr. Then, by the assumption 02.811 and Remark 12.21 it is 
directly shown that 


A ^(c ‘^f){x) >0 for a; G 11. 

Using 03.25p . one can easily verify that 

[ A~^{c~‘^f){x) ■ ^o{x - y) dy > 0, x e OBr, 


'Br 


which together with 03.23p immediately implies that 


(3.25) 

(3.26) 


c = c. 


(3.27) 
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Finally, by fl3.27p and fl3.2ip . one has 

f = I 

The proof is complete. 


□ 


Remark 3.1. We note that the assumption fl2.8p is only needed to guar¬ 
antee that (cf. fl3.24p and fl3.26p i 


A (c f){x) ■ - y) dy x e Br. 


(3.28) 


’Br 


By using a completely similar argument to the proof of Theorem 12.31 
one can show that if / is nonpositive on but strictly negative on 
an open subset. Then fl3.28p also holds true. But as remarked in 
Remark 12.21 fl2.8p is a practical condition. The condition fl 2 . 10 p is 
more general, and this can be observed by simply taking v = \. 


Proof of Theorem 2.4 The proof follows from a similar argument to 
that of Theorem 12.31 Without loss of generality, we can assume that 
/ and Ch are independent of the variable x^, and write them as Cb{x') 
and f{x') for x' = {xi,X 2 ) G By contradiction, we assume that 
there exists another TAT/PAT conhguration {f{x'),c{x')) supported 
in with c{x') given as 

T-\x^) = cf{x^)+T\i:, (3.29) 

where 7 is a positive constant, such that 

Af^c{x,t) = Aj~{x,t) for (x, f) G dVL x M+. (3.30) 

Noting that c~^f — c~^f is independent of the variable 0 : 3 , we have 
by Theorem 12.21 that 

c-V = c-^l (3.31) 

Similar to fl3.24p . we have from fl2.6p in Theorem 12.11 that 


/ ( 1 -c ( 2 /))-A (c f){y)-^o{x-y)dy 

JBr 

= [ (1 - ■ ^~\c-^f){y) ■ *^>o(a; - y) dy 

JBr 


(3.32) 


for any x G Br. Using a completely similar argument in proving fl2.5p 
in the proof of Theorem 12.11 one can show from fl3.32p that 


{l-c-\y))-A-\c-^f){y)-v{y)dy 


'Br 


(3.33) 


= / {l-d-\y))-A-\T-^f){y)-v{y)dy 


'Bf, 


for any harmonic function v. 
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Next, by combining fl2.9p . fl3.29p . fl3.3ip and fl3.33p . we have 


(7 ^ — 7 y A ^(c '^f){x)-v{x)dx = 0 (3.34) 


for any harmonic function v(x). Noting fl2.1UI) . we immediately have 
from (^M> that 


which in turn, together with fl3.3ip . implies that 


/ = /• 


The proof is complete. 


□ 
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